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SYMBOLS  AND  NOTATIONS 
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N 


n. 

1 


r . 
1 


g 


signal  to  be  detected 
Fonrier  transform  of  s 

(2M+l)xl  noise  vector  (white,  Gaussian  and  nncorrelated 
between  any  two  elements) 

the  ith  element  of  N 

received  signal  at  the  ith  sensor 

beam  ontpnt 


D 


variance  of  the  random  quantity  z 
distance  between  two  sensors 


0 


a.  s„ 

c 

A+ 

Y+  (=A+Z) 

L  (={1..}) 


incidence  angle  of  the  incoming  signal  wave  to  the  linear 
array  of  sensors 

sin0,  sin0, 

signal  propagation  speed  in  the  medium 
generalized  inverse  of  the  matrix  A 

minimum  norm  least  squares  solution  of  the  set  of  linear 
equations  AY=Z 

matrix  of  the  low-pass  filter  operator  with  cutoff 

frequency  f  T: 

c 


i  sin[2n( i-j ) fcT] 

1ij  =  jr(i-j)  ;  i»j=0.±l,±2,  ... 

note  that  its  spectrum  is  one  in  the  frequency  range 
(~f  T.f  T) .  and  zero  elsewhere 


( 2M+1 ) x«  matrix  of  the  truncation  operator 

p.  ,  /  li 
ij  l  0;  otherwise 

sampling  interval 


weighting  function 
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CHAPTER  1 
INTRODUCTION 


w 


_ >In  the  area  of  signal  processing,  considerable  attention  has  been 

devoted  to  digital  array  processing  in  recent  years.  This  attention  is 
dne  to  the  increasingly  vide  nse  of  array  processing  for  both  civilian 
and  military  purposes.  Digital  beaaforming,  for  example,  is  popular 
because  of  its  advantages  in  speed,  accuracy,  etc.,  over  conventional 
analog  beaaforming.  Many  papers  concerning  digital  beamforming  have 


recently  been  published  [l]-[8],  It  has  been  shows-  [1],  [8]  that  the 
quality  of  performance  of  a  beaaformer,  such  as  bean  pattern,  signal-to- 
noise  ratio,  etc.,  depends  to  a  great  extent  on  the  number  of  sensors 
used,  i.e.,  the  more  sensors  used,  the  better  the  bean  pattern  becomes. 
In  a  practical  situation,  however,  the  number  of  sensors  say  be 
restricted  by  economical  reasons  or  physical  restrictions .  In  this 
situation  one  may  weight  the  output  of  each  sensor  before  beamforming. 
This  helps  to  some  extent,  although  the  improvement  is  rather  limited. 

A  totally  different  issue,  signal  extrapolation,  has  also  been 
drawing  a  great  amount  of  interest  recently,  largely  in  the  area  of 
spectral  estimation.  It  has  been  shown  that  a  known  portion  of  a  signal 

can  be  extrapolated  outside_of  the  observation  interval  if  the  signal 

\ 

possesses  certain  property  [9],  [12]'.  Many  algorithms ,  both  iterative 
and  non-iterative,  have  been  proposed  for  both  continuous  and  discrete 
cases.-  References  [9]-[13]  serve  as  a  good  review .on  this  issue. 


r  -- 
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The  purpose  of  this  thesis  is  to  use  spatial  signal  extrapolation  in 


digital  beamforming  to  improve  the  beam  pattern  without  adding  more 

O — '  Q 

physical  sensors.  Effectively,  the  sensors  are  added  "’’synthetically” 
through  signal  processing.  This  has  the  potential  to  improve  performance 
considerably.  Other  currently  used  techniques  such  as  interpolation, 
weighting,  etc.  can  be  combined  with  extrapolation  to  achieve  overall 
better  performance. 

\ 
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CHAPTER  2 
BACKGROUND 

Tie  following  discussion  assumes  a  uniform  linear  array 'as  shown  in 
Fig.  2.1.  For  convenience ,  we  assume  that  there  is  an.  odd  number  of 
sensors  in  the  array.  These  sensors  are  indexed  fran  -M  to  H  so  that  the 
total  number  is  2M+1.  Also  we  assume  that  the  signal  is  band-limited  and 
the  noise  is  white  and  Gaussian.  Furthermore,  we  assame  that  the  beam 
steering  specifications  are  met  either  by  sampling  the  sensor  outputs  at 
a  sufficiently  high  frequency,  or  by  using  digital  interpolation 
beamforming  techniques,  as  described  in  [2],  These  assumptions  will  be 
used  throughout  this  thesis. 

2.1  Beamforming 

The  task  of  detecting  a  signal  and  determining  its  direction  can  be 
accomplished  by  digital  beamforming,  i.e.,  by  delaying  and  snming  the 
corresponding  sensor  signals.  Suppose  the  signal  is  a  plane  wave  which 
impinges  upon  the  linear  array  of  2M+1  sensors  at  am  angle  80  with  speed 
c,  (see  Fig.  2.1).  The  received  signal  at  the  ith  semsor  is  denoted  by 

iDa 

rt(t)  -  s(t  +  — -}  (2.1) 

co 

where  a0=Sin90,  The  beamformer  samples  the  output  of  each  sensor,  delays 
the  samples  by  delay  intervals  determined  by  the  steering  angle,  and  then 
sums  these  delayed  samples  of  the  sensor  signals  along  the  array 
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dimension.  This  process  is  described  by 


1  M 

g(a,mA)  -  -  l  ri(mA  - 

2M+1  tt_M  1V  c  1 


■  "r~r  l  .(»A-io(i-!£)] 
2M+1  i— M  c  co 


rr  f  S(e3“"^)  [— i-  I 

**  -it  2M+1 


-i  /  s(.^)yW£--»)].3" 

— Tf  C  C0 


jcjrnA 


din 


(2.2) 


where  Sle^1*®^)  is  the  discrete  Fourier  transform  of  s(mA),  and 
called  the  beam  pattern,  is  given  by 


W(v) 


_ L_  l  e~jivD 

2M+1  i— M 

.  (2M+1)  _ 

sin  vD 

■■■■■■■■  ■'  ■■  ■■■"  e 

(2M+1)  ain  ~ 

2 


(2.3) 


The  magnitude  of  a  beam  pattern  is  usually  displayed  on  a  dB  scale 
as  shown  in  Fig.  2.2  where  the  horizontal  axis  is  i/D.  Notice  that  the 
beam  pattern  is  periodic  with  a  period  of  2n.  It  is  easy  to  see  from 
Eqn.  (2.3)  that  the  largest  value  occurs  at  vD=0 ,  which  coincides  with 
when  the  beam  is  steered  at  the  angle  of  the  incident  signal 


wave . 
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\  v-  v*.  ‘‘T  •" 


Tie  beam  output  at  this  angle  is  g(a0 ,mA)=s(mA) .  Tins  tie  signal  and  its 
incidence  angle  can  both  be  detected.  Notice  tbat  Da/c  must  be  an 
integer  multiple  o£  A  because  of  tie  sampling,  although  Da9/Co  i* 
continuous.  This  implies  that  a  digital  beamformer  cannot  steer 
continuously  in  the  detecting  space.  However,  as  we  mentioned  earlier, 
it  is  assumed  that  the  sampling  rate  is  high  enough  so  that  all  required 
steering  angles  are  covered. 


Also,  it  is  seen  that  the  first  zero  of  W(y)  occurs  when  . 
That  is,  the  main  lobe  width  is  inversely  proportional  to  the  total 
number  of  sensors  in  the  linear  array.  By  increasing  the  number  of 
sensors,  the  width  of  the  main  lobe  of  the  beam  pattern  is  narrowed,  thus 
improving  an  important  measurement  of  the  beam  pattern. 

A  weighting  sequence  w(i)  can  be  applied  to  each  sensor  output 
before  summation.  In  this  case,  Eqn.  (2.2)  takes  the  form 


g(a,mA) 


I  w(i)  rj_[mA 


(2.4) 


and  the  beam  pattern  becomes 


W(v)  =  — -  l  w(i)  e 


2M+1  i=— M 


1  n 

-  [v(0)  +2  1  w(i)  cos  ivDj 

2M+1  i=l 


(2.5) 
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if  v(i)  is  symmetric.  The  weights  w(i)'s  sre  chosen  to  improve  various 
parameters  of  the  ontpnt  beam,  in  ranch  the  same  way  an  FIR  filter  is 
designed  using  a  window  function  to  shape  its  spectrum.  In  general* 
windowing  is  a  well  developed  subject  in  the  literature  on  both  filter 
and  antenna  design  [14],  [15]. 

In  the  case  when  there  is  a  noisy  background,  the  signal-to-noise 
ratio  (SNR)  at  the  beam  output  serves  as  an  essential  measurement  of 
performance.  Assuming  additive  white  Gaussian  noise  which  is 
uncorrelated  from  one  sensor  to  another,  the  received  signal  at  the  ith 
sensor  is 

iDa 

ri<c)  -  s(t  + - -J  +  n^t),  (2.6) 

co 

and 


the  beamformer  output  becomes 


At  the  angle  i-il, 
c  c# 


g(aA) 


s(mA)  y 
21+1  i~M 


v(i)  + 


1 

2M+1 


M  iDa 

I  v(l)  ni(mA - 2.). 

i— M  co 


(2.8) 


The  signal-to-noise  zstio  of  g(mA)  is  defined  as 


M  M  iDa 

SNR  -  10  log  {var[i^l  £  w(i)]  /  var[— -  £  w(i)  ^(mA- 2-}  ] } 

ZM+1  i>i-K  2M+1 


10 


M  .  o  M  , 

'  [  w2(i)} 


logio  K  bb  £  w(i)J  /  z 

2M+1  i— M  (2M+1)2  i— M 


1°  login  -4+10  log  {[  £  w(i)]2  /  £  v2(i)}. 


’10  2 
c 
n 


(2.9) 


i=-M 


The  first  term  in  the  last  expression  is  the  input  SNR.  The  second  term 
is  called  "array  gain".  It  equals  101og10 (2M+1)  when  all  weights  equal 
one.  In  other  words,  beamforming  increases  SNR  by  the  amount  of  array 
gain.  Also,  note  that  the  array  gain  increases  as  the  number  of  sensors 
increases . 


2.2  Extrapolation 

The  motivation  to  link  digital  beamforming  and  signal  extrapolation 
is  quite  obvious  from  the  above  discussion.  Performance  of  the  beam 
pattern  and  beam  SNR  are  dependent  upon  the  number  of  the  sensors  used. 
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Often,  including  a  large  number  of  aensora  in  an  array  ia  not  practical 
becauae  of  the  physical,  geoaetrical,  or  econoaical  reaaons.  An 
excellent  approach  to  solving  thia  problea  ia  veil  known  in  the  area  of 
SAB  (aynthetic  aperture  radar),  where  a  Bowing  aensor  (a  aicrowawe  radar 
carried  on  an  airplane  or  a  satellite)  is  used  to  fora  a  long  synthetic 
array.  Thia  approach  is  not  very  practical  in  sonar  and  underwater  sound 
signal  processing  for  aany  reasons,  such  as  low  propagation  speed,  low 
vehicle  speed,  and  instabilities  in  traveling  along  a  straight  path  [17]. 
Thus,  it  is  natural  to  looh  for  alternate  techniques  which  extend  the 
array  beyond  its  actual  physical  length. 

In  1975,  Fapoulis  [9]  proposed  an  iterative  procedure  for 
extrapolating  an  (continuous)  analytic  signal  based  on  the  observation  of 
only  a  tine  limited  portion  of  it.  Later,  a  series  of  papers  concerning 
this  problea  was  published  ([10]-[13],  etc.).  An  important  paper 
discussing  discrete  extrapolation  is  due  to  Jain  et  al.  [12].  For  the 
purpose  of  this  study,  only  Jain’s  one  step  approach  will  be  discussed 
here . 

Let  the  infinite-dimensional  vector  T={yk;  denote  a 
sequence,  and  the  (2M+l)xl  vector  Z  denote  an  observation  of  T  over  a 

limited  interval,  i.e.,  2j*yj  for  Using  the  truncation  operator 
P,  we  have 


Z=PT . 


(2.10) 


If  yk'$  are  the  samples  of  a  continuous  function  y(t)  which  is  band- 
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V 


(Note  that  the  subscripts  here  are  different  from  the  conventional  row- 
colnmn  index  of  a  matrix),  and  HI  is  a  sqnare  matrix  of  dimension  (2M+1) 
with  entries  given  by 
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B 

K 

i> 

i*', 

i 

I 


tin[2»(i-j)fcT] 

hlij  ■  — <2-15> 

la  cate  of  noisy  observation 

2-PY+N,  (2.16) 

ve  have  a  "mean-square  extrapolation": 

I«H(Hl+yI)“*Z  (2.17) 

where  y*«r*/<y*  and  1  is  the  identity  matrix  (also  see  [12]). 

i  A 

The  extrapolated  sequence  Y  (also  Y)  does  not  equal  Y,  in  general. 
However,  it  converges  to  Y  (hence  to  y(t))  when  the  interval  T  approaches 
zero  (One  more  condition  is  needed  for  this  conclusion,  i.e.,  y(t)  has 
finite  energy,  see  [13]}. 


I 


/ 

/ 


CHAPTER  3 


EXTRAPOLATION  BEAMFORMING 

3.1  Real-Time  Implementation 

Jain's  one-step  extrapolation  procedure  can  be  used  for  spatial 

extrapolation,  as  veil  as  time  extrapolation.  Specifically,  suppose  the 

x-axis  describes  spatial  position  along  tbe  linear  array,  and  a  plane 

wave  s(t+x  — )  impinges  upon  this  array  witb  an  incidence  angle  0a 
co 

and  speed  c0  as  previously  shown  in  Fig.  2.1.  The  received  signal  at  the 
ith  sensor  is  s(t+iDjp)  as  in  Eqn.  (2.1).  Let  the  signal  s(t)  be  band- 
limited  over  (“fc,fc).  Hence  if  S(f)  is  the  continuous  Fourier  transform 
of  s(t),  then  S(f)“0  for  |fl>fc.  For  a  fixed  t,  the  traveling  wave 
s(t+x|^)  is  a  function  of  x.  Without  loss  of  generality,  we  may  as  well 
assume  that  t=0 .  Then  we  define 


u(x)=s(^x)  . 
co 


(3.1) 


Its  Fourier  transform  is  given  by 


U(k)=|2S(|5k) .  (3.2) 

*0  *0 

Since  S(j^k)=0  for  we  find  that  u(x)  is  spatially  band-limited 

& 

over  (~c^fc , c^fc) .  Furthermore,  by  the  well  known  sampling  theorem,  in 
order  to  reconstruct  u(x)  from  its  samples  u(iD),  the  interval  D  must 
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D  < 


(3.3) 


where  i*  the  minimum  wave  length  of  the  signal.  Since  a«ssind«  is  no 
greater  than  unity,  we  get  fired  bounds 


(3.4) 


and 


D<^.  (3.5) 

2 


Then  u(r)  is  (spatially)  sampled  at  a  rate  1/D,  the  sequence  u(iD)  is 

band-limited  over  (-1  D,k,D) ,  or  (-— f„ ,^-f.) .  Thus,  Jain's  approach  can 

c  c  c  0  c  C  q  w 

be  applied  spatially  to  each  group  of  samples  obtained  at  2M+1  sensors  at 

each  temporal  saapling  time.  The  T's  in  Eqn.  (2.14)  and  Eqn.  (2.15)  are 
D 

replaced  by  — 's  with  D  satisfying  the  same  requirement  as  conventional 

beamformer,  i.e.,  D<  — .  A  real-time  implementation  scheme  is  shown  in 

2 

Fig.  3.1.  The  truncation  point  M'  of  the  extrapolation  along  z-azis  will 
be  discussed  in  section  3.3.  From  the  figure  we  see  that  in  eztending  the 
array  length  synthetically,  we  add  some  multipliers,  adders.  and 
registers.  This  is  worthwhile  even  in  comparison  with  the  case  where 
physical  sensors  are  added.  This  is  because  with  the  compatible  results 
(discussed  later),  the  synthetic  approach  avoids  possible  long  noise- 
sensitive  connection  cables  between  the  additional  sensors  and  the 
central  processing  unit,  and  saves  sensors  and  A/D  converters  as  well. 


/ 
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Other  techniques,  such  ss  temporal  interpolation,  snap  ling 
techniques  for  band-limited  signals,  and  frequency  domain  beaaforaing 
[2]-[7],  can  be  nsed  simultaneously  with  the  spatial  extrapolation. 
However,  they  will  not  be  discussed  in  this  thesis. 

3.2  Extrapolation  Error 

As  mentioned  earlier,  the  extrapolated  sequences  do  not  in  general 
equal  the  samples  of  the  original  continuous  signal  which  we  seek  to 
reconstruct,  i.e.,  there  is  some  "extrapolation  error".  Te  also 
mentioned  that  these  extrapolated  sequences  converge  to  those  continuous 

values  uniformly  if  the  sampling  interval  T  (—  in  our  case)  approaches 

ce 

zero  [13].  Thus,  for  a  relatively  small  —  we  may  expect  that  the 

c» 

extrapolation  error  is  within  some  bound  so  that  the  extrapolation 
beaaforming  is  meaningful.  The  quantitative  measurement  of  the  error  is 
needed  in  evaluating  the  performance  of  the  extrapolation  beanformer. 
Unfortunately,  it  is  very  difficult  to  evaluate  these  errors  analytically 
for  the  reason  discussed  below. 

Ve  notice  that  Jain's  formula  was  obtained  earlier  through  numerical 
approximation  of  some  integral  equations  by  Cadzov  [10].  Specifically, 
the  one-step  extrapolation  for  continuous  case  is 

•9- 

f  hl(c-T)  x(t)  dT  -  z( c) ;  t«A  (3.6) 

T€A 

y(t)  -  /  h(t-T)  x(t)  dt;  t^A 
T€  A 


(3.7) 
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where  h(t)  end  hl(t)  ere  the  inpulse  response  of  e  low-pass  filter  with 
cutoff  frequency  f  x(t)  is  an  intermediate  variable,  and  A  is  the 
observation  interval.  One  first  solves  Eqn.  (3.6)  for  x(t)  based  on  the 
observation  z(t),  then  substitutes  the  x(t)  in  Eqn.  (3.7)  to  obtain  the 
result  y(t).  A  numerical  approximation  to  this  set  of  equations  is 

M 

l  hl(lT-JT)  x( jT)  -  z(iT);  |  i  |  <  M  (3.8) 

j— M 

M 

y(lT)  **  l  h(lT-jT)  x(jT)  T;  |i|>M.  (3.9) 

j— M 

Note  that  this  set  of  equations  is  exactly  the  same  as  Eqn.  (2.13)  where 
it  was  written  in  matrix  form.  Since  Eqn.  (3.6)  and  Eqn.  (3.7) 
reconstruct  the  continuous  signal  exactly,  the  extrapolation  error 
mentioned  above  is  just  the  numerical  approximation  error  of  Eqn.  (3.6) 
and  Eqn.  (3.7).  We  now  show  that  it  is  impractical  to  evaluate  this 
error.  First  of  all,  Eqn.  (3.6)  may  not  have  a  solution  at  all  if  y(t) 
(hence  z(t))  does  not  satisfy  certain  conditions  [10],  [11].  In  this 

case,  although  an  HNXS  solution  to  Eqn.  (3.8)  exists,  it  is  meaningless 
to  mention  the  error.  These  conditions  will  not  be  quoted  here,  since 
even  if  a  solution  to  Eqn.  (3.6)  does  exist,  we  still  cannot  estimate  the 
error.  This  is  because  the  numerical  approximation  to  Eqn.  (3.6) 
(solution  of  Eqn.  (3.8))  can  be  quite  different  from  the  samples  of  the 
continuous  solution  of  Eqn.  (3.6)  due  to  the  ill— posed  nature  of  the 
Fredholm  equation  of  the  first  kind.  Equation  (3.6)  can  be  rewritten  as 
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M 

/  hl(lT-T)  x(t)  dT  -  T  l  hl(iT-JT)  x(jT)  +  Q(hl,x,T) 
t«A  j— M 


-  z(lT);  |  1  |  <  M 


(3.10) 


where  Q(hl,x,T)  denotes  the  non-zero  reaainder.  Equation  (3.8),  however, 
is  the  above  equation  with  Q(hl,x,T)*0.  Thus,  Equ.  (3.8)  actually  solves 
another  integral  equation 

M 

/  hl(il-T)  x(t)  dt  -  T  l  hl(iT-jT)  x( jT)  +  Q(hl,3C,T) 
t<tA  j— M 

-  z(iT)  +  Q(hl,S,T);  |i|<M  (3.11) 

with  Q(hl,x,T) *0 .  By  the  argument  in  [16],  we  see  that  the  small 
perturbation  Q(hl,x,T)  on  the  right  hand  side  of  Eqn.  (3.11)  causes  the 
solution  x(r)  to  move  arbitrarily  far  from  desired  x(x) .  Hence,  the 
error  x(jT)-x(jT)  cannot  be  predicted,  let  alone  the  error  of  y(jT)  in 
Eqn.  (3.9) . 

A  practical  approach  to  this  problem  is  to  conduct  sotie  numerical 
experiments  with  example  signal  sequences,  and  then  to  model  the  error 
obtained  from  the  experiments.  This  has  in  fact  been  done,  and  the  main 
results  are  presented  below.  The  method  used  to  invert  the  matrix 
(HI+7I)  in  Jain's  formula  is  the  well-known  Gauss  elimination  method 
which  gives  more  accurate  and  more  reliable  results  than  other  methods 
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such  as  Levinson's  algorithm.  The  computation  time  is  aot  important 

bscsass  tbs  coefficients  of  H(m+yl)“x  need  to  be  cslcalsted  only  once 

before  the  installation  of  the  equipment.  A  sum  of  two  shifted  sine 

functions  (called  SHSINC)  and  a  sine  function  (SINE)  are  used  as  signals. 

Although  the  latter  has  not  been  proven  to  converge  to  continuons 

solutions,  Jain's  approach  still  gives  a  good  results  as  shown  below. 

The  observation  is  made  at  17  points,  i.e.,  M**8.  The  quantity  f  ~  i* 

c  c  ^ 

£ 

the  digital  cutoff  freqnency  F,  which  is  less  than  0.5  (since  D< j^~) • 

c 

The  digital  signal  frequency  f— ,  denoted  by  FI  (for  SHSINC,  it  is  the 

ce 

highest  frequency  of  the  sine  function), is  less  than  or  eqnal  to  F. 

First,  we  investigate  the  noise-free  case,  i.e.,  N-0.  In  a 

practical  situation,  the  signal  sequences  may  have  several  degrees  of 
freedom.  Their  frequencies  as  well  as  their  phases  may  vary.  We  now 
look  at  the  extrapolation  error  for  different  phase  shifts.  It  is 
observed  for  SINE  shifting  from  0  to  2ji.  For  SHSINC,  the  separation 
points  between  two  main  peaks  is  60,  the  shift  is  from  the  -150th  point 
to  the  50th  point.  The  filter  cutoff  frequency  is  0.05  if  it  is  not 
specified.  Dashed  lines  are  signals,  and  solid  lines  are  error  if  not 
otherwise  indicated. 

Figure  3.2  and  Fig.  3.3  show  the  extrapolation  error  vs.  x/D  for 
SHSINC  and  SINE  respectively  when  their  shifts  are  all  zero.  We  can  see 
a  flat  region  in  the  central  part  of  the  extrapolation  error,  indicating 
that  good  extrapolation  almost  doubles  the  observation  length.  Notice 
that  7  is  not  zero  as  it  should  be  in  ideal  noise-free  case.  This  is 
because  experimentally,  very  poor  extrapolation  results  if  7=0,  due  to 
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Fig.  3.2  Extrapolation  error  vs.  i  for  SHSINC 
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Fig.  3.3  Extrapolation  error  vs.  i  for  SINE 


the  ill-conditionness  of  the  matrix  El  and  the  round  off  error  in 

computation.  This  problem  is  solved  by  adding  a  very  small  number  to  the 
diagonal  of  El  to  "stabilize"  it  [12],  which  is  exactly  the  same  as 
letting  y  equal  this  small  number  (in  our  case  10“1X  to  10~X*  is  enough) 
in  Eqn.  (2.17)  for  N=0. 

Next  we  fix  i  (i-x/D)  at  a  point,  let  the  ‘  signal  shift  in 

aforementioned  ranges,  and  observe  the  extrapolation  error  vs.  shift  for 
this  particular  point  of  i.  The  plots  for  SHSINC  are  shown  in  Fig.  3.4 
through  Fig.  3.7  for  different  i  and  different  FI.  It  is  interesting  to 
see  that  the  error  is  almost  a  shifted  version  of  the  signal  with  the 
magnitude  multiplied  by  a  negative  number.  It  is  also  observed  that  this 
model  fits  better  for  low  signal  -frequency  case  (Fl=0.02).  Figure  3.8 
through  Fig.  3.13  show  same  plots  for  SINE.  We  observe  that  in  these 
cases  the  error  is  exactly  sinusoidal  with  one  complete  cycle.  As  FI 
increases  to  F,  its  magnitude  becomes  larger,  and  its  shift  with  respect 

to  the  signal  becomes  irregular  (see  Fig.  3.12  and  Fig.  3.13  ,  the  shifts 

0 

are  not  180  as  it  is  for  smaller  FI  case)  . 

The  magnitude  of  the  peak  error  for  these  shift  ranges  at  each  point 

of  i  is  shown  in  Fig.  3.14  through  Fig.  3.18  for  different  cases  (since 

this  magnitude  is  symmetric,  only  one  half  is  shown).  It  is  clearly  seen 
that 

(1)  For  the  signal  of  SINE,  the  error  becomes  very  large  when  FI  equals 
F  while  it  remains  small  even  for  F1=0.96F.  This  phenomenon  is  not 
observed  for  SHSINC.  This  is  because  the  spectrum  of  SINE  is  an 
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Fig .  3.4 


SHSINC 
C  2M  -  1  )  - 
I  -50 
F-  .  05 
F  1  -  .  0  4 
SNR  -110 


17 


Extrapolation  error  vs.  shift  for  SHSINC,  y-10~XX ,  Fl=0.04,  x*30D 


Fig.  3.5  Extrapolation  error  vs.  shift  for  SHSINC,  -y=10-il 
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Fig.  3.6  Extrapolation 
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error  vs.  shift  for  SHSINC,  y=*10-11,  Fl=0.02,  x=30D 
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Fig.  3.7  Extrapolation  error  vs.  shift  for  SHSINC,  -y=10 
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Fig.  3.8  Extrapolation  error  vs.  shift  for  SINE,  y=10-11,  Fl=0.02,  x-30D 
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Fig.  3.9  Extrapolation  error  vs.  shift  for  SINE,  y=10-li,  Fl=0.02,  x=50D 
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Fig.  3.10  Extrapolation  error  vs.  shift  for  SINE,  7=10-11,  Fl=0.04,  x=30D 


Fig.  3.11  Extrapolation  error  vs.  shift  for  SINE,  y=10  11,  Fl=0.04,  x=50D 
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Fig.  3.13  Extrapolation  error  vs.  shift  for  SINE,  ^=10  lt,  Fl=0.05,  x=50D 
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Fig.  3.14  Magnitude  of  extrapolation  error  ▼*.  i  for  SHSINC,  y=10  ,  F=0.05 
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Fig.  3.15  Magnitude  of  extrapolation  error  vs.  i  for  SHSINC,  y=10 
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iapulse,  and  if  it  equals  the  cutoff  frequency  of  the  low-past 
filter,  it  is  alaost  cut  off  by  the  filter.  While  the  apectrua  of 
SHSINC  extends  from  zero  to  FI,  thus  no  inforaatiou  is  lost  even  if 
FI  equals  F. 

(2)  For  a  fixed  FI,  saaller  filter  frequency  F  (but  still  larger  than 
FI)  achieves  better  extrapolation  in  the  sense  that  the  error  is 
saaller  (e.g.,  compare  the  curve  of  Fl*0.04  in  Fig.  3. Id  with  the 
curve  of  same  FI  in  Fig.  3.17). 

(3)  The  extrapolation  error  tends  to  unity  when  i  increases  as  shown  in 
Fig.  3.18.  This  is  because  the  extrapolated  sequence  tends  to  zero 
as  i  increases. 

Since  SINE  signal  is  a  basis  for  any  other  signals  via  Fourier 
transfora,  and  also  exhibits  a  better  property  in  extrapolation  in  the 
sense  that  the  error  vs.  shift  is  an  exact  sine  function,  we  will 
concentrate  on  the  extrapolation  and  beaaforaing  problems  for  the  sine 
signal  only. 

Suppose  the  sine  signal  has  a  random  phase  9  which  is  uniformly 
distributed  on  [0,2jt],  i.e.,  y(x)*sin(Qx+0) ,  then  by  the  above 
observation  we  can  model  the  error  at  a  fixed  i  as 

E(x,0)=E( iD,0)=-A( i) sin(2iD^9)  (3.12) 

so  that  the  error  and  the  signal  are  180  out  of  phase  for  each  x.  The 
magnitude  function  A(i)  can  be  obtained  by  curve-fitting  the  plots  of 
magnitude  vs.  i.  It  is  easy  to  calculate  the  density  function  of  E  at 
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From  observation  (1)  above  and  from  the  earlier  observation  abont 
error  vs.  shift  we  see  that  the  filter  cutoff  frequency  f  (hence  F) 
should  be  larger  than  maximum  possible  signal  frequency  f  (FI).  On  the 
other  hand,  observation  (2)  istplies  that  it  is  better  that  F  is  small. 
Thus,  in  a  practical  situation,  we  need  to  estimate  the  maximum  signal 

frequency  f  and  select  ffi  tsd  D  (F=fc— )  carefully  to  trade  off  between 
these  two  aspects. 

Now  let  us  consider  the  noisy  case.  For  uncorrelated  white  Gaussian 
noise  N  added  to  the  observation,  we  have 

z=Zi+N  (3.14) 

where  Z1=py.  and 

I  =  H(Hl-HrI)"1Z 
*  HdLL+yl)  (Z1+N) 

=  H(Hl+TrI)~1Z1+H(Hl-HrI)_1N.  (3.15) 

Thus  T  can  be  expressed  as  a  sum  of  two  terms 
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$■ 


V*x 


(3.16) 


where 


^■H(H1+y1)_1Zi  (3.17) 

^-HfHl+YD^N.  (3.18) 

A  A  A 

Notice  that  Y#  Y^  Y^  are  all  inf inite^diaenslonal  rector*  rhile  Zi 

A 

tad  N  are  (2M+l)xl  vectors.  Ve  see  that  each  element  of  Y^  is  a  linear 

A 

combinatioa  of  that  of  N,  heace  Yr  i#  also  white  aad  paussian.  The 
analysis  for  is  exactly  the  same  as  before  except  that  y  is  ouch 

larger.  An  arbitrary  value  for  y  is  chosea  so  that  101ogloi  is  13  dB 
(This  valae  will  be  used  throughout  the  thesis  if  it  is  aot  otherwise 
specified) . 

The  plots  for  error  vs.  phase  shift  are  shown  in  Fig.  3.19  through 

o 

Fig.  3.28.  Ve  see  that  the  error  coincides  with  the  signal  (with  180 
/ 

out  of  phase)  better  than  N=0  case  for  FI  less  than  0.8F.  Beyond  F1=F, 
the  error  does  aot  fit  into  the  model  well.  The  error  magnitude  vs.  i  is 
plotted  ia  Fig.  3.29  and  Fig.  3.30  for  0^i^50  (A(i)  is  symmetric).  We 
see  also  that  the  error  gets  larger  for  FI  beyond  0.8F.  Thus,  in 
designing  an  extrapolation  beamforaer,  we  should  choose  ffi  such  that 


J.8f 

c 


>f 


max  • 


(3.19) 
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Fig.  3.23  Extrapolation  error  vs.  shift  for  SINE,  y*0.05,  Fl=0.02,  x*15D 


Fig.  3.24  Extrapolation  error  vs.  shift  for  SINE,  7=0.05,  Fl=0.02,  x-25D 
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Fig.  3.25  Extrapolation  error  vs.  shift  for  SINE,  y-0.05,  Fl-0.02,  x*50D 
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Fig.  3.26  Extrapolation  error  vs.  shift  for  SINE,  y=0.05,  Fl=0.04. 
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Fig.  3.27  Extrapolation  error  vs.  shift  for  SINE,  f=0.05,  Fl=0.048,  x»21D 
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Fig.  3.28  Extrapolation  error  vs.  shift  for  SINE.  y=0.05.  Fl=0.05,  x=20D 
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Thai  is  these  two  plots,  discarding  the  carves  for  F1)0.8F.  the  remaining 
three  corves  (F1*0,  F1=*0.4F  and  F1=0.8F)  are  averaged  and  are  plotted 

with  the  dashed  lines.  They  will  be  osed  is  Section  3.3  and  is  Chapter  4 
as  the  magnitude  function  A(i)  in  the  error  model  (Eqn.  (3.12)). 

Notice  that  in  Fig.  3.1*  the  inpnt  to  the  beamformer  for  lil{X  is 
fed  directly  from  the  sensors  without  extrapolation.  This  is  because  we 
want  to  use  as  much  real  information  as  possible.  Thus  the  elements  of 
for  I i |<M  are  replaced  by  those  of  Z  and  A(i)  is  set  to  zero  for  ii](M. 
This  will  not  affect  the  extrapolation  for  |i|>M. 


3 .3  Evaluation  Of  Perforaanc< 


At  this  stage  we  are  ready  to  evaluate  the  performance  of  the 
digital  extrapolation  beamformer  using  the  error  model  obtained  in  last 
section.  Specifically,  for  a  sinusoidal  incoming  wave 


»(t,x)  =■  sin[uQ  (t  + — -)  +  $] 


(3.20) 


with  a  random  phase  uniformly  distributed  over  [0,2;;],  the  ith  sensor 


receives 


ri(t,iD)  =  sin[u0  (t  + - — ]  +  $]  +  n^t); 


|  i  |  <  M 


(3.21) 


where  n^t)  is  the  aforementioned  noise.  After  temporal  stapling,  tb.es 

V*  fora  the  observation  vector  Z,  or  the  elements  of  f  (denoted  yp  f0 
lila.  i.e.. 


?i(mA,iD) 


iDa 

sin[oi0(mA  + - -)  +  $]  +  n^mA); 

co 


j i  |  <  M.  (3.22] 


Then  they  are  extrapolated  to  |i|=M'#  i.e.. 


iDa 


yi(mA,iD)  =  sin[a)0(mA  + -—■)  +  *]  -  A(i)  sin[cuQ  (nA  + - 2.]  +  *] 

+  ^(mA); 

M  +  1  <  U  |  <  M'  (3.23) 


where  ^(mA)  is  the  extrapolated  noise,  i.e.,  the  ith  element  of  the 
vector  for  M+lN<IiI^M'. 

All  these  2M'+1  elements  of  Y  are  fed  to  the  beamformer  which  gives 
an  output  as  follows: 


g(a,mA)  =  - 

2M’+1 


a 

l  sin{w0[nA-iD  f— - 2.)]  +  *1 

i 31  “M  *  c  c0  1 

iv» 

[  A(i)  sin{u  [aA-iD  f— - } 

i=i(M+l)  c  c0" 


l 

2M'+1 


wit*  re 
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^  f  n^mA);  j  i  |  <  M 

UiCmA)  -  |  ftl(aA).  Mfl  <  f  1  J  <  M* 


(3.25) 


and 


i-±(M+l) 


M’  -M' 

I  +  I 

i-M+1  i— (M+l) 


(3.26) 


We  now  ewalnate  two  major  measurements :  beam  pattern  and  signal-to- 
noise  ratio. 

3.3.1  Beam  Pattern 

Witbont  considering  tie  noise,  Eqn.  (3.24)  becomes 


.  M’  a 

g(a,mA)  -  — — —  I  sin{u0[raA-iD  (— - -)  ]  + 

ZM  +1  -  - 


c  c. 


±M’ 


~  ~T77~T  l  AC1)  sin{u0[niA-iD  (— - -j  ]  +  t>] 

"  i*±(M+l)  c  c0 


M' 


_A  f  s(ej«mAj  r  1  £  fi-  jwiD(i  -  +  j<?  — 


2tt  1 


2M'+1  i=-M' 

±M' 

r 

l 

A(i)  e^101^ 

=(M+1) 

dw 


(3.27) 
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where  Sfe^*1*^)  ie  the  diecrete  Fourier  transform  af  sinw,aA _  Tith  V 
denoting  «(§-|^),  the  negnitnde  of  the  bees  pattern  is 


jw(v)  |  - 


1 

2M’+1 

1 

2M’+1 


I  «'13DV  -  I  A(i) 
i— M’  i-±(M+l) 


M'  M' 

1  +  2  J  cos  iDv  -  2  £  A(i)  cos  iDv 

i-1  i-M+1 


(3.28) 


where  the  symmetry  property  of  A(i)  is  used.  Ve  also  get  a  closed  form 
for  the  first  summation,  i.e. , 


|W(V) 


1 

2M’+1 


1_e-j(2M»+l)Dv 


ejM'Dv 


tM’ 

l  A(i)  e"1JDU 
i— —  (X+l) 


2M  '+1 


sin  IS™?.  Dv 


sin 


Dv 


M* 

“2  £  A(i)  cos  iDv 

i-M+1 


(3.29) 


A  particular  case  for  Eqn.  (3.29)  is  when  M'-M  (no  extrapolation),  and 
Eqn.  (3.29)  degenerates  to  Eqn.  (2.3).  When  MOW,  the  beam  pattern  is 

modified  by  a  factor  -  2  ^  A(i)  cos  IDv  besides  the  expected 

i-M+1 

improvement  (M  being  increased  to  M').  This  is  of  coarse  caused  by  the 
extrapolation  error,  and  results  in  some  degradation  of  the  beam  pattern. 
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Vhen  a  window  (weighting)  is  added  in  beaaforaing,  the  beam  output 
(Eqn.  (3.24))  becomes 


,  M'  a 

g(a,mA)  -  — — — -  l  v(i)  ain[oi0(oA-iD(~  ~  -—■) )  +  <fr] 
2M'+1  ±m-H>  c  co 


2M'+1  i-±(M+l) 
M* 


M*  a 

l  w(i)  A(i)  sia[(0o(mA-iD(-2-  -  — — •) )  +  <fr] 

c  c0 


.  1 


2M'+1 


I  w(i)  'Sji  (mA  - 
i-M'  c 


(3.30) 


where  the  window  w(i)  is  2M'+1  point  long,  i.e.,  and  is 

symmetric.  The  corresponding  beam  pattern  magnitude  is 


|w(v) 


1 

2M*+1 


L  w(i)  e  J 
i— M  ’ 


-  I  w(i)  A(i)  e  1^°y 

i-±(M+l) 


1 

2M’+1 


M’  M* 

w(0)  +2  £  w(i)  cos  iDv  -  2  £  w(i)  A(i)  cos  iDv 

i-1  i-M-rl 


(3-31) 


Also,  it  degenerates  to  Eqn.  (2.5)  when  M'=M  and  degenerates 
Eqn.  (3.28)  when  w(i)=l  for  each  i,  a  rectangular  window. 


to 
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3.3.2 


Let  ns  reconsider  Eqn.  (3.25).  The  n^ 'g  are  nncorrelated  for  each  i 
**  assumed.  Also,  the  ft^'s  are  some  linear  combination  of  n^'a.  Let  hjj 
denote  the  element  of  the  (2M'+l)x(2M+l)  matrix  Hdl-tyl)""1,  at  the  ith  row 
and  the  jth  column  where  -M'<i<M',  Thus  n^(mA)  can  be  expressed 


*»i(mA)  -  l  j  nj(mA);  JH-1  <  |  i  {  <  M* . 


(3.32) 


Then  n^nA)  can  be  expressed  in  terms  of  n^(mA) : 


ni(mA); 


I  i  I  <  M 


"i(mA)  -  •  M  (3.33) 

1  hjj  nj(mA);  M+l  <  |  1  J  <  M». 
j— M 


Notice  that  the  n\(mA) 's  are  no  longer  nacorrelated  for  M+14  I  i  |(M' , 
At  the  searching  angle  when  |»|*.  Eqn.  (3.24)  becomes 


1  r 

g(a0,mA)  -  sin(ui0nA-h*i)  -  — - -  sin(u)0nA-t~Ji)  T  A(i) 

2M'+1  i-t(M+l) 


2M'+1 


Mf  iDa 

v  ~  r  .  Oi 

L  n^raA  -  — — ■; 


(3.3A) 
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The  variance  of  j(»0,mA)  is  • 


^MV 

var[g(a0,mA)]  -  c2  [l - - —  £  A(i)] 
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(3.35) 


where  aj  ij  the  variance  of  the  signal,  the  sane  as  <j*  on  page  11. 
The  signs 1-to-noise  ratio  at  the  bean  output  is  then 


SNK  -  10  log10  {o2s  [l  - 


1 


2M'+1 


±M'  2 
l  A(i) ]  }  -  10  log, 

i=~(M+l) 


10 


(2M’+1) 
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+M’  M  ±M*  ±M’  M 

+  2  l  I  £  I  I^kJ  ^j]} 
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M 
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(3.36) 


where  the  symmetry  property  of  A(i)  is  again  used.  In  the  particular 
case  when  M'»M,  all  the  summation  terms  no  longer  exist,  and  Eqn. 
(3.36)  degenerates  to 


2 

a  , 

SNR  -  10  logl0  -§“  +  10  logig  [  (2M+1)  ]  -  10  log1Q  (2M+1) 
an 


2 

a 

-  10  log,  Q  -§-  +  10  log10  (2M+1) ,  (3.37) 

a 


which  coincides  with  Eqn.  (2.9)  for  a  rectangular  window. 


For  M’>M,  we 
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M  ,  ±M» 

[  I  »J(D  +  2  I 


M 
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(3.39) 


Tien  the  signal-to-aoise  ratio  at  the  beaa  output  is 
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M 
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(3.40) 


Note  that  this  expression  reduces  to  Eqn.  (3.36)  for  a  rectangular  window 


and  reduces  to  Eon.  (2.9)  for  M'=M  as  expected. 


CHAPTER  4 


EXPERIMENTAL  RESULTS 

In  this  chapter  the  results  of  digital  computer  simulation  will  be 
presented.  Some  terminology  used  in  the  figures  and  in  the  discussion 
for  this  chapter  is  listed  below.  The  observation  points  (number  of 
sensors)  are  17  throughout  this  discussion,  i.e.,  M«8,  as  in  Section  3.2. 
Also,  the  filter  cutoff  frequency  F  will  be  0.05  if  it  is  not  otherwise 
specified. 

Input  SNR 
Beam  SNR 
' xdB  BW 
xdB  SMER 

SLML  R 

Peak  Loss 

4.1  Main  Results 

Figure  4.1  shows  a  plot  of  Eqn .  (3.29),  beam  output  SNR  vs.  M’,  the 
extrapolation  length,  for  two  cases  when  the  input  SNR  is  30  dB  and  13 
dB,  respectively.  At  the  beginning  point  of  these  curves  where  M'=M=8 
(no  extrapolation),  the  beam  output  SNR  is  42.3  dB  and  25.3  dB  for  each 


signal-to-noise  ratio  at  each  sensor,  dB 
signal-to-noise  ratio  at  beam  output,  dB 
xdB  beam  width  of  the  main  lobe 

side  lobe  -  main  lobe  energy  ratio  when  those  above 
xdB  are  taken  as  in  main  lobe  and  those  below  are 
in  side  lobes,  dB 

side  lobe  -  main  lobe  magnitude  ratio,  where  the 
largest  side  lobe  is  used,  dB 

main  lobe  peak  loss,  dB 


i 
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Fig.  4.1  Bean  SNR  vs.  extrapolation  length,  F»0.05,  7-0.05 
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Fig.  4.2  Beam  SNR  vs.  extrapolation  length,  vith 
Hamming  windows,  F=0.05,  y=0.05 
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cafe.  In  other  words,  beamforming  increases  SNR  by  12.3  dB.  It  is 
interesting  to  see  that  the  beaa  SNR  displays  some  decaying  oscillation 
as  M'  increases  and  tends  to  some  certain  value.  At  tbe  beaa  SNR 
achieves  a  local  maximum  which  is  41.2  dB  for  30  dB  input  SNR  and  24.3  dB 
for  13  dB  input  SNR.  These  values  are  only  one  dB  less  than  those  at 
M**8.  At  the  double  length  of  the  array.  JT»17,  the  beaa  SNR  gets  the 
miniauo  value  39.0  dB  and  22.0  dB.  respectively.  Figure  4.2  shows  the 
same  plot  when  Hamming  windows  are  applied  for  each  M* .  In  this  case  the 
oscillation  is  smoothed  out.  The  array  gain  at  M'«8  becomes  10.8  dB. 
The  maximum  value  occurs  at  M'*15  where  the  values  happen  to  be  42.3  dB 
and  25.3  dB.  the  same  as  those  at  M'»8  in  Fig.  4.1.  And  the  values  at 
M'=29  equal  their  starting  values  40.8  dB  and  23.8  dB,  respectively. 
From  these  plots  we  see  that  the  beam  SNR  is  not  improved  by 
extrapolation,  in  fact  it  decreases  slightly.  This  is  probably  due  to 
the  strong  correlation  of  the  n^'s  and  the  n^'s. 

We  now  observe  the  beam  patterns.  Figure  4.3  shows  the  beam  pattern 
for  M ' =M=8 ,  no  extrapolation  case  for  rectangular  windows.  The  plots  for 
M'=17  and  H'=29  when  extrapolation  is  applied  are  shown  in  Fig.  4.4  and 
Fig.  4.5.  We  observe  that  the  extrapolation  narrows  the  main  lobe  width 
extensively.  Moreover,  it  decreases  the  side  lobe  -  main  lobe  magnitude 
ratio  and  the  side  lobe  -  main  lobe  energy  ratio.  The  peak  loss  is  not 
that  important  since  the  beam  SNR  decreases  little  for  these  M"s  as 
shown  in  Fig.  4.1.  Of  course  M'=29  is  preferred  since  the  performance  of 
the  beam  pattern  is  slightly  better  than  that  in  M'=17  and  the 
degradation  of  SNR  is  only  one  dB.  As  a  contrast,  the  beam  pattern  for 
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foil  observation,  l.t..  oil  the  input*  to  the  beaaforaer  ore  obtoiaed 
froa  season  (M'-M),  is  shown  ia  Fig.  4.6  for  M-29.  It  is  seea  that  the 
besa  width  ia  this  esse  is  even  Borrower  then  that  in  Fig.  4.4  os  we 
expected,  siace  aoae  of  these  59  inputs  to  the  beoaforaer  contains 
extrspolotion  error.  However,  the  side  lobe  -  a* in  lobe  aagnitude  ratio 
and  the  side  lobe  -  a* in  lobe  energy  ratio  reaaia  the  soae.  Thus,  in 
trading  off  aaong  these  aeasnreaents,  the  digital  extrapolation 
beaaforaer  is  coapatible  with  a  conventional  digital  beaaforaer  which 
nses  double,  even  triple  the  number  of  sensors.  It  is  also  observed  that 
the  iaproveaent  in  the  performance  does  not  increase  rapidly  as  M' 
increases.  In  fact  beyond  about  M’-40,  the  side  lobe  -  main  lobe  ratio 
starts  to  increase  (not  shown  here) .  Thus,  in  consideration  of 
complexity  of  implementation  and  the  iaproveaent  of  the  performance, 
M'*29  is  approximately  optimum  in  this  case. 

When  Hamming  windows  are  applied  in  the  digital  extrapolation 
beaaforaer,  results  are  obtained  as  shown  in  Fig.  4.7  through  Fig.  4.10. 
Figure  4.7  shows  the  case  when  no  extrapolation  is  applied  to  the  17 
sensors.  Comparing  it  with  Fig.  4.3,  it  is  clear  that  the  application  of 
Hamming  windows  decreases  the  side  lobe  -  main  lobe  ratio  on  the  expense 
of  widening  the  main  lobe.  The  results  of  applying  extrapolation  are 
still  good  as  shown  in  Fig.  4.8  and  Fig.  4.9.  The  main  lobe  is  narrowed 
while  the  side  lobe  -  main  lobe  magnitude  ratio  is  further  reduced  than 
the  cases  without  applying  Hamming  windows.  The  extrapolation  length 
M'=15  is  preferred  because  of  its  smallest  side  lobe  -  main  lobe 
magnitude  ratio  (-27  dB)  and  largest  beam  SNB  (see  Fig.  4.2).  For  M'=29 
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(Fig.  4.9).  the  aide  lobe  -  main  lobe  magnitude  ratio  ia  only  -23  dB» 
thus  ve  would  bare  used  tbe  eaae  shown  ia  Fig.  4.5  where  much  aarzower 
aaia  lobe  overwhelms  tbe  compatible  aide  lobe  -  aaia  lobe  magnitude  ratio 
with  tbe  former.  Figure  4.10  aerrea  agaia  aa  a  coatraat  wbieb  ia 
obtaiaed  aaiag  fall  observation  of  31  aeaaora.  We  aee  that  $he  aide  lobe 
redactioa  for  tbe  eztrapolatioa  beaaforaer  ia  aot  as  good  as  that  for  tbe 
conveational  beaaforaer  ia  applying  Hamming  windows. 

4.2  Effects  Of  Paraeter  Changes 

As  we  aentioned  earlier,  tbe  paraaeters  F  and  y  can  be  changed  to 
alter  tbe  performance.  To  gain  aore  insight  into  this,  we  now  observe 
tbe  case  when  F-0.1.  Figure  4.11  and  Fig.  4.12  show  the  beaa  SNR  for 
rectangular  windows  and  Hamming  windows  as  before.  They  show  same 
character  as  that  in  Fig.  4.1  and  Fig.  4.2  except  that  tbe  oscillation 
frequency  is  increased  ia  Fig.  4.11.  Two  peaks  nearest  to  tbe  origin 
along  M'-axis  are  at  M’~19  and  M'»29.  Tbe  beam  patterns  at  these  points 
are  shown  in  Fig.  4.13  and  Fig.  4.14  (The  beam  pattern  for  M ' =M=8  is  same 
as  that  in  Fig.  4.3  since  no  extrapolation  is  involved  when  M'=M) .  We 
see  that  all  measurements  are  not  as  good  as  that  shown  in  Fig.  4.4  and 
Fig.  4.5  while  they  are  still  much  better  than  those  in  Fig.  4.3.  The 
case  when  Hamming  windows  are  applied  is  shown  in  Fig.  4.15  and 
Fig.  4.16  for  M’=13  (where  beam  SNR  gets  the  largest  value)  and  M'=19  (in 
contrast  with  Fig.  4.13).  The  same  observation  is  made  as  the  case  with 
the  rectangular  windows.  It  is  also  observed  for  both  cases  that  the 
improvement  in  the  beam  pattern  increases  more  slowly  as  M'  increases. 
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Fig.  4.13  Bean  pattern  for  17  sensors,  extrapolated  to  39,  F-0.1,  y»0.05 
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Fig.  4.14  Beam  pattern  for  17  sensors,  extrapolated  to  59,  F=0.1,  y=0.05 
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Fig.  4.15  Beam  pattern  for  17  sensors,  extrapolated 

to  27,  with  a  Hamming  window,  F«0.1,  y«0.05 


BEAM  PATTERN  CDS) 


( 2M-  !  )  -  1 7 
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1  209  BW-  . 5027 
3C3  SMER--S.S34 
1  0CB  S  E R  -  -  1  4 . 7i 
SLMLR- -21 .02 
-EA<  L0SS-7.53P 


Fig.  4.16  Beam  pattern  for  17  sensors,  extrapolated 

to  39,  with  a  Hamming  window,  F=0.1,  y=Q.0S 
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*•  conclude  from  thut  observations  that  smaller  extrapolation  filter 
frequency  F  performs  better.  Tbie  coincide*  with  tbe  criterion  (2) 
(Observation  (2))  in  Section  3.2.  Uni  it  is  important  to  select  D  and 

*c  ‘“el1  that  Eqn.  (3.3),  Eqa.  (3.19)  and  tbe  criterion  (2)  in 

Sectioa  3.2  are  aet  in  designing  a  digital  extrapolation  beamformer. 

All  tbe  experiments  shown  above  were  conducted  under  tbe  assumption 

that  tbe  parameter  y  in  Eqn.  (2.17)  is  chosen  such  that  101og1#I  i*  3.3  dB 

o*  * 

even  if  tbe  input  SNR  (-lOlog^-f)  i,  aot  13  dB.  This  is  because,  for 

XL 

one  reason,  tbe  hardware  implementation  (see  Fig.  3.2)  requires  that  y  be 
fixed.  Besides,  tbe  input  SNR  is  usually  not  known  in  a  practical 
litutioni  this  also  requires  a  good  guess  for  y.  Tbe  value  for  y  that 
has  been  used  can  then  be  changed.  A  group  of  plots  is  shown  in 

•  4.17  through  Fig.  4.20  for  beam  SNR  and  beam  pattern  as  we  did 
before,  where  r  is  chosen  such  that  101ogl#I  is  30  dB.  Ve  observe  the 
similar  behavior  with  the  previous  ones  while  some  differences  are  also 
noticed.  Tbe  beam  SNR  performs  worse  in  the  sense  that  it  tends  to  some 
lower  value  (even  lower  that  input  SNR)  and  decreases  very  fast  as  M' 
departs  from  8.  On  the  other  hand,  however,  the  beam  pattern  performs 
better  in  the  sense  that  the  main  lobe  width  is  narrower  and  the  side 
lobe  -  main  lobe  ratio  is  smaller  compared  with  the  previous 
corresponding  ones.  Ibis  gives  ns  another  degree  of  freedom  in  choosing 
parameters  to  compensate  among  various  design  requirement  and 
measurements  of  performance. 


various 


BEAM  SNR  (OB) 
50  r 


.  INPUT  SNR 

-  50  DB 

.  INPUT  SNR 

-  15  DB 


Fig.  4.17  Beam  SNR  vs.  extrapolation  length.  F-0.05,  7=0.001 


WITH  HflM.WNDWS 
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-  50  DB 

2 .  INPUT  SNR 
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Fig.  4.18  Beam  SNR  vs.  extrapolation  length,  with 
Hamming  windows,  F=0.05,  7=0.001 
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Fig.  4.19  B«ta  pattern  for  17  sensors,  extrapolated  to  35,  F-0.05,  7>-0.001 


BEAM  PATTERN  (OB) 


f  2M- 1  )  - ! 7 
( 2M '  -  1  )  -35 
WITH  HAM.WNOW 
30B  BW- . 267 
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I 0DB  SMER- -  15. S 
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Fig.  4.20  Beam  pattern  for  17  sensors,  extrapolated  to  35, 
with  a  Hamming  window,  F=0.05,  y=0.001 
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CHAPTER  3 
SUMMARY 

la  tbit  thesis  the  application  of  spatial  extrapolation  in  digital 
beamforming  was  studied.  An  essential  problea  in  evaluating  performance 
of  the  digital  extrapolation  beamformer  is  to  model  the  extrapolation 
error.  Section  3.2  was  devoted  to  this  issue  for  both  finite  energy 
signals  and  periodic  signals.  An  error  model  was  obtained  as  a  random 
function  with  certain  distribution  for  a  particularly  important  kind  of 
signals,  sinusoidal  input  signals.  Having  accomplished  this,  it  is  not 
difficult  to  construct  a  real-time  hardware  implementation  scheme. 

A  large  amount  of  digital  computer  simulation  was  done  for  various 
cases.  The  results  are  quite  promising.  It  was  shown  that  the  beam 
pattern  can  in  general  be  improved  extensively,  in  some  cases  we  can  even 
extrapolate  as  far  as  the  double  length  or  triple  length  of  the  physical 
sensors  with  the  performance  compatible  to  a  conventional  digital 
beamformer  with  this  large  amount  of  sensors.  On  the  other  hand, 
however,  signal-to-noise  ratio  is  degraded  slightly  by  the  extrapolation. 
A  parameter  y  can  be  chosen  to  compensate  between  these  two  measurements. 

Other  techniques  such  as  windowing,  temporal  interpolation,  etc.  can 
be  combined  with  the  spatial  extrapolation  to  improve  performance 
further.  Examples  using  Hamming  windows  were  included. 

Since  the  sinusoidal  wave  is  a  basis  function  for  other  signals 
through  Fourier  series  and  Fourier  integral,  and  the  one-step 
extrapolation  process  is  linear,  the  results  obtained  in  this  thesis  can 


/ 
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be  readily  extended  to  tbe  general  case.  On  the  other  hand,  however, 
sinnsoidal  wave  contains  infinite  energy.  The  convergence  of  its 
discrete  extrapolation  to  its  continuous  value  has  not  been  proven  yet 
[13].  This,  together  with  a  theoretical  sudel  for  its  extrapolation 


error,  remains  open  questions  in  the  digital  extrapolation  beaaforaing. 
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